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1.  Introduction 

We  consider  in  this  paper  a two-dimensional  diffusion  process  that 
arises  in  conjunction  with  the  following  type  of  aueuing  system.  The  system 
is  composed  of  two  single  server  facilities  (or  stations)  arranged  in  tandem. 
Customers  arrive  individually  from  outside  the  system  and  uueue  up  I'or 
service  at  the  first  station.  Having  completed  service  there,  they  pi-oceed 
to  a queue  in  front  of  the  second  station,  and  after  completing  service 
at  the  second  station  they  depart  the  system.  Bervice  is  by  order  of 
arrival  at  each  station. 

We  assume  that  the  system  is  empty  at  time  " and  that  the  inter- 
arrival times  and  service  times  at  the  two  station^  „ .i  three  mutually 
independent  seciuences  of  IID  random  variables.  The  interarrival  times  are 
distributed  as  a random  variable  v^  having  mean  a^  and  finite  variance 
Vq,  and  the  service  times  at  station  k (k  l.P)  are  distributed  as  a 

random  variable  v,  with  mean  a,  and  finite  variance  | '.  . We  define 
k k k 

the  traffic  intensities  (k  ■ 1,P).  The  system  is  said  to 

be  stable  if  < 1 for  each  station  k,  and  a situation  of  heavy 
traffic  is  said  to  prevail  if  . ^ is  near  one  for  each  station.  F'or  a 
stable  system  we  further  liefine 


r 


(( 


(a^^-a^l  ^ A (i->. 


Thus  a stable  system  is  in  iieavy  traffic  when  u is  positive  but  near 
zero.  Finally,  it  will  be  convenient  to  define 


^'I 

‘^1  1, 

M ■ 

and  - 

bp 

where  = 1q  Ty  = - r^,  o',,  --  ^ 

- (a^-a^)/a.  Note  that  >1^  and  t are  the  mean  vector-  and  covariatu-e 

matrix  respectively  of  the  random  vector  Because  we  do 

9- 

not  require  that  the  variances  be  strictly  positive,  .|.  neeil  not  he 

positive  definite,  but  it  is  non-nej^.ative  definite  in  any  case. 

Let  be  tiie  waiting;  time  (exclusive  of  service  time\  at 

n 

station  k for  the  n^*^  arriving  customer,  and  let  W = (W^,  W . If 

n n n 

the  system  is  stable,  it  is  well  kno-.vn  that  there  exists  a random  vector 
1 

W - (W  ,vr  ) such  that  W ==>  W as  n -><«,  where  -->  lienotes  weak 

n 

convergence  (convergence  in  distribution). 

In  an  earlier  paper  (2l,  a limit  theorem  was  proved  to  show  that 


r ’ 

P P 2 

the  case  Yq  = = Y'o)  invariance  principle  argument,  but  the 

general  problem  of  evaluating  the  limit  distribution  was  left  open.  It 
is  to  this  problem  that  we  now  return.  The  current  discussion  is  restricted 
to  the  case  of  two  queues  in  tandem,  but  most  of  what  will  be  said  extends 
readily  to  the  general  model  with  K stations  in  series. 

In  Sections  2 and  5 we  repeat  the  definition  of  Z and  show  that 
it  is  a diffusion  process  (continuous  strong  Markov  process)  whose  state 
space  is  the  non-negative  quadrant.  More  specifically,  on  the  interior 
of  its  state  space  the  process  behaves  like  ordinary  Brownian  Motion 
with  drift  vector  ^ and  covariance  matrix  and  it  reflects  instan- 
taneously at  each  axis  (boundary  surface).  At  one  axis  (the  one  corresponding 
to  an  empty  second  station  in  the  original  nueueing  system) , the  direction 

t 

! 

j of  reflection  is  normal  to  the  axis,  but  at  the  other  axis  the  reflection 

; has  a tangential  component  as  well.  In  Section  ^ we  also  display  whe 

' (weak  infinitesimal)  generator  of  Z. 

r 

In  Sections  through  6 we  demonstrate  that  the  limit  distribution 
of  Z is  the  solution  of  a first  passage  problem  for  a certain  dual 
diffusion  process  Z*.  Using  the  analytical  theory  of  Markov  processes, 

I 

[ we  then  derive  in  Section  7 the  following  partial  differential  equation 

I (with  boundary  conditions)  for  the  density  f of  the  limit  distribution. 

' ^1^11  '^12^12  '*■  I °2^22  " ^•‘l^l  ■ " 0 , 

'°12^L7  ^ } 4'^,'  " M;f](x,0)  - 0 , 
f|  a^f^  + (|  t-  - ^aj_f](0,y)=  0 . 

L - 


(Here,  and  Later  in  tlie  paper,  we  use  the  notation 


fl(x,y)  - ^ fU,y) 


and  so  forth.'  In  Section  H we  show  that  the  limit  distribution  has 
independent  components  iff  vdiic.h  case  each  component  is 

exponentially  distributed.  This  expands  slightly  tiie  class  of  solutions 
obtained  previously  in  [I’J.  Finally,  in  Section  9,  we  explicitly  solve 


for  f when  is  the  identity  matrix  and 


0,  (in  tenr.is  of  the 


original  tanden  cueuing  system,  this  corresponds  to  the  case  = p.^, 
deterministic  services  at  the  first  station,  and  = P.-,-  ■ The  distri- 
bution doer,  not  have  independent  components. 

The  explicit  solutions  obtained  here  are  certainly  f.''agmentary . 
but  we  hope  that  they  can  help  to  stimulate  and  direct  the  interest  of 
others  in  this  problem.  To  determine  e.xplicitly  the  general  solution 
of  the  partial  differential  equation  displayed  above,  one  must  overcome 
some  significant  analytical  difficulties.  Relatively  little  effort  has 
been  devoted  to  the  problem  thus  far,  however,  and  an  explicit  general 
solution  may  in  fact  be  obtainable. 

Finally,  we  note  that  tandem  queues  are  the  simplest  case  of  a 
queuing  network,  so  any  analytical  insights  obtained  here  should  be 
valuable  in  determining  heavy  traffic  approximations  for  the  general  net- 
work model.  Limit  theorems  .justifying  diffusion  approximations  for  general 
networks  will  be  presented  in  a forthcoming  paper  ny  Reiman  |1)], 
together  with  analytical  characterizations  of  the  relevant  diffusions. 


2.  A Reflection  Mapping  in  Two-Dlmenslons 


Let  B denote  the  set  of  functions  x(-)  • fx^(-).  x (•)/ 
that  map  10,«)  into  B and  are  bounded  on  finite  intervals.  Let 

be  the  subset  of  such  functions  that  are  non-negative  and  non-decreasirif^ 
For  X B let  C(x)  be  the  set  of  all  t C B^  such  that 

Xj^(t)  + _0  for  all  t ^ 0 

and 

x^(t)  - ^^(t)  + P.^lt)  > 0 for  all  t _ 0 . 

It  is  obvious  that  C(x)  is  non-empty.  Moreover,  it  has  a minimal  element, 
meaning  that  there  exists  y ..  C(x'  such  that  1 'j(t.  a.nd 

yp(t)  < ?,,(t)  for  all  f,  C(x)  and  all  t > 0.  This  minimal  element  is 

explicitly  given  by 

y (t)  = sup  (-X  (u)]’^ 

0 < u < t 

(1) 

yp(t)  = sup  {-[x^,(u)  -y^(u)]]'^=  sup  { I -x  ivi  l"*"-  x,^(u)1^  . 
0<u£t  “ 0<v<u_t  ^ 

We  define  a mapping  f ^ (f^,fp);B  B^  by  setting  y^  f^(x  anJ 
y,,  - f,(x).  Also,  we  define  g = (g.|^,g,  i:B  by  setting  gJ^x^=x^^y^ 

and  g,,(x)  ^ Xp  - y^  + y,,. 

An  interp>-etation  of  the  mappings  f and  g is  as  follows. 

Imagine  a man  (hereafter  called  the  controller)  who  monitors  two  bank 
accounts,  the  contents  of  which  fluctuate  due  to  uncontrollable  deposits 


and  withdrawals.  In  the  absence  any  intervention  by  the  controller, 
the  content  of  account  k at  time  t will  be  • 1“  particular, 

the  initial  contents  of  the  two  accounts  are  x^(0)  and  x^(0).  The 
controller  has  the  ability  to  transfer  money  from  account  ■ to  account  1, 
and  we  denote  by  §^(t)  tlie  total  amount  of  money  so  transferred  during 
the  interval  [0,t].  He  also  can  deposit  money  (gotten  outside  the 
system)  in  account  2,  and  we  denote  by  ^^(t)  the  total  amount  deposited 
during  [0,tj.  The  pair  i = is  called  a control.  If  we  suppose 

that  the  controller  is  required  to  keep  the  content  of  each  account  non- 
negative, then  C(x)  represents  the  set  of  feasible  controls.  If  we 
suppose  that  the  controller  wants  to  minimize  the  total  volume  of  his 
transactions  (in  the  strongest  possible  sense),  then  y = f(x)  represents 
the  optimal  control,  and  the  content  of  account  k at  time  t will  be 

z (t)  when  the  optimal  control  is  employed,  where  z - g(x). 

K 

The  effect  of  the  reflection  mapiping  g is  porti'ayed  graphically 
in  Figure  L below.  Suppose  x £ B is  continuous  and  represents  the  sample 
path  of  a stochastic  process.  Let  y = f(x)  and  z = g(x) . Continuing 
the  interpretation  of  f and  g given  above,  one  might  call  x the 
uncontrolled  process  and  z the  controlled  process.  If  x(0)  is  strictly 
positive,  then  x = z until  the  first  time  that  x strikes  a boundary 
point  of  the  non-negative  quadrant.  If  the  horizontal  axis  (x  , 0) 

is  struck  first,  then  the  local  behavior  of  the  controlled  process  after 
hitting  is  given  by  - x^  and  z^  = x..,  + y^.  The  control  increases 

in  the  minimal  amounts  necessary  to  keep  z,,  non-negative,  and  we  have 
the  normal  reflection  pictured  in  Figure  1.  If  the  vertical  axis  (x^  - 0) 


Figure  1 


Angles  of  Reflection  Induced  by  the  Mapping  g 

a aonick  first,  however,  then  the  effect  of  the  mapping  g is  slightly 
more  complicated.  The  local  behavior  of  the  controlled  process  after 
hitting  is  now  given  by  ~ ^ ^2  ~ ^ control 

y^  increases  in  the  minimal  amounts  necessary  to  keep  non-negative 

and  we  have  non-normal  reflection  in  the  direction  shown  in  Figure  1. 

A final  important  property  of  the  mapping  g is  that  it  is 
memoryless  in  the  following  sense.  Suppose  x c.  B anu  z = g(x).  For 
arbitrary  T > 0,  let  z,j,(t)  - z(T  + t)  and 

x^(t)  = z(T)  + fx(T  + t)  - x(T)] 


for  t > 0.  Then 


(2) 


= g(X,j.)(t) 


for  t > 0 . 


r r 


Thus  the  sample  path  of  tlie  controlled  or  reflected  process  alter 

time  T is  completely  detemiined  by  7.(T)  and  the  sajuj)Je  path  of  the 
uncontrolled  process  x after  time  T.  To  be  specific,  it  is  gotten 
by  applying  the  ref'ection  mapping  g to  the  process  X,j,.  The  memory  less 
property  {2)  follows  immediately  from  the  chai-acterization  of  f(x  as 
the  minimal  element  of  C(x)  but  is  not  at  all  appai-ent  from  the  explicit 
representation  (l). 


t>.  The  Diffusion 

Viewing  ^ and  ^ as  arbitrary  data  hereafter  (except  for  the 
requirement  that  be  symmetric  and  non-negative  definite),  let 

X = (X(t),  t > 0)  be  a two-dimensional  Brownian  Motion  with  drift  vector 

, '^2 
covariance  matrix  and  general  starting  state.  If  o'^  *0  ^'l' 

2 2 2 2 2 

cr,,,  = -Y-.  and  cj„  = (see  section  I),  then  such  a i)rocess  can 

JJr  X £Z  X 

be  constructed  by  taking 

t > 0 , 

where  the  processes  ^ are  independent  standard  (zero  mean  and  unit 

K 

variance!  Brownian  Motions  v;ith  general  starting  states.  In  tlie  usual 
way,  we  denote  by  ^y('^  distribution  on  the  path  space  of  X 

corresponding  to  starting  state  x - (x^,x  ) C ]R  ‘ . We  define  the 


■ X^(t) 

Xp(t) 

r p,(t)  - + X t 

-Li  c:  (L  i 

H 


processes  Y - f'.X)  anO  '/  k(X).  so  that, 

Z^(t1  - X^(ti  + Y^(t)  . t _ 0 , 

Z,(t)  = X^(t)  - Y,(t'  ^ Y^(t)  , t ^0  . 

Proposition  1.  Z is  a strong;  KarXov  process  with  state  spat  e S IR4  , 
and  its  sample  paths  are  almost  surely  continuous. 

Proof.  From  the  explioih  representation  for  y - f(^x)  that  wa.'^  given 
in  Section  2 it  follows  easily  that  f and  g both  map  continuous 
functions  into  continuous  functions.  Thus  the  path  continuity  of  7. 
follows  from  the  continuity  of  X.  The  strong  Markov  property  is  immediate 
from  (2)  and  the  strong  Markov  property  and  stationary,  independent  incre- 
ments of 

Observe  that  Z(0)  = (O)  whenever  X(0)  . S.  Thus 

can  alternately  be  interpreted  as  the  distribution  on  tlie  pati:  space  of 
Z corresponding  to  ZfO)-  x,  provided  that  x G 3.  In  tlie  interest  of 
characterizing  more  completely  our  underlying  diffusion  Z,  we  conclude 
this  section  with  a result  shov.'ing  how  its  generator  operates  on  a class 
of  smooth  functions.  We  shall  make  no  further  use  of  the  generator  of 
Z in  this  paper,  however. 

Following  I\ynkin  |1J,  v?e  define  ^ to  be  the  set  of  bovinded. 
measurable  functions  f; S ->  .IH  such  that  T^f  f as  t i 0.  wliere 
T.f(x)  - E [f(Z(t))|.  We  denote  by  & the  set  of  f <5^  such  that 

L>  X 

(T^f  - f)/t  converges  boundelly  pointwise  as  t I 0 to  a limit  in  . 


The  limit  is  denoted  Gf,  the  operator  G thus  defined  is  called  the 
weak  infinitesimal  (generator  of  Z ''or  just  genei'aior  foj-  .dioi’t  ■ , 
and  b'  is  called  the  domain  of  the  generator. 

Let  5 denote  the  set  of  bounded,  continuous  functions  1!< 

such  that  all  partial  derivatives  of  f exist  and  are  boon  le  1 and 
continuous  on  S.  Let  (for  smooth)  be  the  set  of  functions  f 
that  vanish  on  some  neighbohood  of  the  origin.  In  saying  mat  f ha.s 
continuous  partial  derivatives  on  the  boundary  of  S,  v/e  mean  that  each 
partial  derivative  approaches  a finite  limit  at  the  boundary  and  the  limit 
is  a continuous  function  of  the  boundai-y  point.  Symbols  like  f.,(xj^,0i 
are  understood  throughout  to  be  shorthand  notations  for  such  limits. 

From  the  path  continuity  of  Z and  bounded  convergence  it  follows 
that  every  bounded,  continuous  f:S  IF  Is  in  ^ . Thu.s 
For  ease  of  notation,  we  define  a differential  operator  P by 


The  following  proposition  will  not  be  proved,  since  it  will  be  subsumed 
by  a more  general  result  of  Reiman’s  [^]  and  it  is  not  used  in  the 
rem.ainder  of  this  paper.  A very  similar  result,  Proposition  < of 
Section  will  be  proved  in  detail,  however. 

Proposition  2.  Suppose  f jJ  . Then  f ^ O'  iff 


(^) 

f^(x^,0)  = 0 

for 

x^  > 0, 

and 

(h) 

(f^-fp)(0,x.,)  - 0 

for 

X,,  > 0, 

in  which  case 

Gf  = 'f  . J . 

10 


Conditions  {^)  nnd  {h)  require  that  at  each  boundary  surface 
(axis)  the  directional  derivative  in  the  direction  siiovm  in  Fii^ure  1 
be  zero.  In  the  definitional  system  of  Watanabe  7|,  such  restrictions 
on  the  domain  of  the  generator  are  used  to  define  direction  of  refiection 
at  a boundary  point.  Thus,  from  the  perspective  of  Watanabe' s theory, 
Proposition  2 justifies  our  statement  that  the  diffusion  Z reflects 
instantaneously  at  its  boundai-y,  the  an(ile  of  reflection  being  as  siiown 
in  Figure  I. 


ij.  The  Limit  Distribution 

Assming  that  the  process  Z begins  at  the  ori^if-.  we  v,lsh  now 
to  characterize  (and  if  possible  compute)  the  disti-lbutiori; 

! F^(z)  = PQ(Z(t)  ^ z)  for  t > 0 and  z S . 

[ 

I 

I It  will  actually  be  rriore  convenient  to  work  wit}i 

= pQ(Z‘]_(t)  < z^,  Z^(t)  + Z^(t'  < t-  z ' . 

r1 

defined  for  t > 0 and  (z,  ,z,d  £ S.  The  reason  is  ■c'ne  followir;g  result, 

— id 

' which  is  precisely  analogous  to  Lemir;a  1 of  ' 

Proposition  j.  i>^(zj^,  -•  P^IM,  (t)  ^ z^,  M,j(t)  < t z^|  . wiiere 

M (t)  = sup  (X  (u))  , 

0<u^t 

M,j(t)  - sup  (X^, (v)  t Xj^(ui)  . 

' 0-_v_u_t  ‘ 

11 

^ l-l-j 


Proof.  The  argument  ts  so  similar  to  the  proof  of  hemma  1 in  jP]  that 
we  sliall  mereiy  sketcli  it.  First  use  fl)  to  write  riown  an  explicit 
representation  for  Z(t)  = (Z^(t),  Z,^(t)),  observing  that  Uie  genera] 
representation  simplifies  considerably  when  X(0)  - 0.  Let  t > 0 be 
fixed.  Assuming  that  X(0)  • 0,  observe  that  {X(u),  0 ^ u < t)  has 

the  same  distribution  as 

{X(t)  - X(t-u),  0 < u < t)  . 

since  X has  stationary  independent  increments.  Maxing  this  substitution 
in  the  representation  foi'  Z(t)  and  simplifying,  we  obtain  the  desired 
proposition. 

With  X(0'  - 0,  it  is  clear  that  and  M,  (•)  are  almost 

surely  non-decreasing,  so  the  limits 


- Lim 
t — 


and 


- dim  M,  (t) 
t 00 


exist  almofjt  surely  Furthermore, 


F^(x)  ->F(x)  and 


as  t , 


for  each  S,  where 


Finally,  I'roin  the  strong  inw  I'or  Brownian  Motion  it  folLov;,:  that 
• ■”  iff  ■;  0,  v;hfjfeas  M , < oj  iff  <'  0 and  + ,,,  <■  0. 

Thus  we  iiave  the  following,. 

Proposition  1.  The  limit  distribution  F is  propei-  iff  , ^ - 0 and 

li-|  '*■  ir,  ■ 

By  definition,  is  tiie  mcucimum  of  the  unj-estrit;1.eii  Brownian 

Motion  X^.  When  |...,  < 0 and  X^(0)  = 0,  is  knovm  to  iiave  an 
exponential  distribution  with  mean  o^/P|(i^|,  and  tliis  gives  tiie  first 
marginal  distribution  of  F.  Tiie  same  result  can  be  obtained  by  obsei-ving 
that  Z,  is  just  modified  by  a reflecting  barrier  at  .:ero.  and  then 

using  well  known  result, s for  ref lecte  i Brownian  Motion.  Note  that  , 
is  not  just  X,  modified  iiy  a reflecting  barrier  at  r.ero.  In  fact,  the 
one-dirnenslonal  process  Z,^  is  not  Markov. 


t.  A Related  Diffu.sion  Prccess 

Let  X'<(t'  ■ - X(t  foi'  t _ 0,  and  let  distri- 

bution on  the  r ath  space  of  X‘  correspenuing  to  X ■ .'0  x for  x F . 
('This  is  tiasier  than  writing';  negative  subscripts.'  Next  let 

Y'(t'i  sup  f -X'(u)  I 
0 < u • t 


and 


infft  _ 0:X’(tl  - Y*(t)  •_  0 ) , 


if  the  indicated  t-set  is  empty.  Now  define  a ];rocess 


with  T*  = oj 


Z'^  = fZ*(t),  t > 0)  by  setting 


Z“(t)  = 


X-(ti  - Y,'^(tJ 


X*(l”'  j - Y^('rA-) 


if  0 < t < T* 


if  t ^ T"  , 


X:<;(t)  + Y'^.ft) 


X*(T')  ^ Y,'<(T^) 


if  0 < t • T>- 


if  t > T' 


Observe  that,  with  Z;  defined  in  terms  of  Y*  in  this  way,  Y*  is  ti»e 

> / ‘ 

minimal  function  of  X**  such  tiiat  Z/(t)  j 0 for  all  t ■ 0.  ilere- 
after  we  restrict  attention  to  paths  of  X><  that  begin  in  the  non-negative 
quadrant,  this  insuring  that  Z''(-)  c.  S. 

The  process  7.-^  is  obtained  from  X'<  in  two  staffer. 

First,  we  impose  an  instantaneously  refiecting  V;arrier  at  the  axis  X'  0, 
the  angle  of  reflection  being  as  shown  in  Figure  . below.  Occond,  an 


Figure 


Angle  of  Reflection  for  the  ] roce. 


absortiin^  barrier  ia  inijioaed  at  the  :ixia  - 0.  The  proof  of  the 
foLLowint’:  is  virtually  identical  to  the  proof  of  !'roi)OS it.ion  1,  so  we 
ielete  it. 


Proposition  'j.  Z*  is  a strong;  Markov  process  wl  tn  .state  .space  15,  and 
its  sample  paths  are  almost  surely  continuous. 


Let  G*  he  the  v;eak  infinitesimal  generator  of  ana  let 

Q'  be  its  domain.  Let  be  ciefined  as  in  Section  j.  Finally,  define 
a differential  operator  D**  by 


D-f  = i c,-' -I 


')  ° -1^10 


-1^1  - 


Proposition  6.  Suppose  f % . Then  f . iff  (f.  ^-f  , Kxj^,0 ; 0 

for  1.-  0.  In  this  case,  G'fi,  /,)  - 0 for  z •-  and  G*f  - D' f 

otherv/ise. 


Proof.  Assume  f L ^ find  let  T|f(.z)  E^l  f (Z'<  (t) ) 1 for  c S and 

t > 0.  We  first  determine  the  conditions  undei'  vdiich  ;T'f  - 1')  ^t 

converges  pointwise  ac  ■‘.10. 

For  each  interior  point  r,  of  S,  let  b^(’.'  demote  the  pi'ob- 

ability  that,  startin,'  in  state  z,  the  boundai’y  i.s  struck  by  Z* 

before  time  t.  Using  nothing  moi-e  td'ian  the  first  exit  disti'ibution 
of  one-aimen:i  icnal  Hrowniar;  Motion  ' i'ecMuse  our  state  space  i.':  a rectangle 
it  is  easy  to  show  tliat  o(t  for  eacti  intiM’ioi  point  z. 


Thu 


T*f(z)  - E,^lf(X«(t))]  + o(t)  , 


because  Z.''(t)  - X’'(t)  until  the  first  hitting  of  the  boundary.  It 
follows  that  |T*f(z)  - f(/.)l  't  ->D'*f(z)  as  t i 0,  since  !)•  is  the 
generator  of  X*  on  (See,  for  example,  section  i.f  of  Ito  lit).) 

If  u (O.z.J  with  z,,  > 0,  then  T^f(z)  f(z)  for  all  t ^ 0. 
since  Z^lt)  = 2 for  all  t >0.  Thus  (T^f(  - f'/.)|  t as  t t 0. 

To  examine  the  boundary  - 0,  we  nee  i some  more  notation.  Tor  purposes 
of  this  proof  only,  let 


M(t)  . sup  lX,^,(u)  - X«(0)  ! , 

0<u<t  ' 

Mx-(t)  • sup  ‘X,^(u)  -X  (Oil  sup  X-(O')  -X>(ui| 
0^u<t  0-_u_t 

for  t ' 0.  Thus  Y<(t)  --  iMxCt;  - X’'(0)  |''’  for  t _ 0.  Ob.sei-ve  that 
the  distributions  of  the  Brovmlan  maxima  M'lt)  and  hM.  10  not  depen  1 
on  X,x(0l.  Furthermore,  the  maximum  distrilmtion  is  kno.-/n  explicitly, 
and  direct  computation  I’eveals  that 

(?)  E*J  M’' (t)  I 't  « and  E,^JM(t)|  t -)■«  f>-s  t.  . 0 . 

(8)  E^[  (jM<(t))'']/t  ->  Up  and  EX;f(M(t))  |/i.  as  t I 0 

Until  further  notice,  let  / (r.^,01  with  2,  ■ u.  It  is  -acU  known 

that  X'/tl  t Mx(tl  then  has.  the  sane  iistrihutiou  a.-:  M(i'  for  eeu.d; 


t > 0.  (See  the  proof  of  Proposition  ^ above.)  Thu.;, 


Ey  (X^(t) t X'(Li  M-(f'  I (M-i't  I E'  ' i , ; . 

Dividing:  both  sides  of  ( i'  by  t,  lettint'’  t i 0.  -tn  1 u.^irif;  ('' 
obtain 

(10)  Eyxy(t'  M"{t)]/t  ->-  - 0 , as  t i 0 . 

Finally,  letting  /X^-'t)  X|^(t)  ~ X<(0;,  we  need 

(11)  E;(h\X^(t'!  M'(t)|/t  as  t i 0 . 

For  the  case  u.  - 0,  this  is  proved  as  foLiow.s.  Fir.-.t  ob: 

(IP)  Ey-'Xl'(t)  M<':t)l  =E»'.X'(t>  (M(t'i  -X'(V('I 

for  t > 0.  (See  again  the  proof  of  Preposition  aljove.  Alsc 
gives  (when  0) 

(ly)  E'i-'xr(t  M-'t)j  ^ Ey/x-(t  M(tii  . 

Z L ‘t’.  I 

Combining  (!,'')  and  [Ij]  gives  (11).  The  {'eneral  cu:;e  i.s  |<r'  vu.i 
it  to  the  driftless  case  through  obvious  bounds. 

As:;uniing  as  above  that  X'.(O’)  - Z,»(0)  - o and  i' (o 
define  'X'(t)  as  before  anl 


, we  thei. 


;erve  that 


, synmietry 


!)y  I'ediioin: 


0 , we 


AZ'(t)  = X'(t)  t M'(t)  . 


(14)  AZ^^(t)  = - M*(t), 

Then 

(1‘,.)  T^f(z)  = E*:  f(zj^  + AZ*(t),  Z/(tV  J + o(t)  . 

the  o(t)  ten;!  on  ttie  right  correspoirting  to  thie  evunt  fT‘  ■ tl.  Now 
an  application  of  Taylor's  Theorem  (with  the  exact  form  of  the  remainder 
gives,  when  combine.]  with  (!'  ' , 

(l6)  T^fiz)  = f('x)  + f^(/)  E<'AZ^(tl  I + ;’.:(t ' ' 

^ I fj^f.)  e;|  (A:-ft))‘  1 + i f ;.-(t.. 

f f^,,(rO  E-’  -Z^vti  . o{t  , 

where  ^ rernnirler  term.  Fnrtt'.ermore,  since  the  Lnii’.i  order 

partial  derivatives  of  f are  bounded,  it  ir  easy  to  that 

H^(z)  = o(t),  using  ,iust  criide  bounds  on  tlie  thii'.i  moments  of  ■.Zj'Ct' 
and  AZ»(t). 

As  v/e  shall  indicate  shortly,  eacli  of  the  ■'•coua-or 'ei'  teimis  on 
the  right  si.le  of  ( 1'  ) converges  to  a finite  Limit  w.heii  '-.e  .livi  ie  it  by 
t and  let  t 1 0.  The  first-order  terms  can  be  n^written.  (using  'll). 
a.s 

f-^(z.)„^t  t-  f (z),.  t ^ If  ,(zi  - 1 E'jM-(t'l  . 


la 


0V  Ti.  II 


Thus  it  follows  from  (?)  that  - f(s)],'t  converges  to  a finite 

limit  iff  f,,('i)  - f^(-'’-)  - 0.  Wlien  this  holds,  we  combine  ( U ) with  (iM 
and  (B)-(ll)  to  obtain 

(in  [T«f(:n  - r(z)]/t  ->D*-f(  -.)  - i ((t'^  + 


as  t i 0.  in  this  calculation,  we  have  also  used  the  o'nvious  relationships 


E^j  (.•  X|(t))^i  .'t 

as  t 4 0, 

Al  (fXJ^(t)  X-'l't)  1 ''t  ->  cr‘^. 

as  t J 0, 

E^l  (X';(t)  )■  1 -t  ct  ' 

as  t i 0. 

Since  (f,  - f ) 0 alori{'  tiie  ;ixis  s = 0,  it  follows  llutt. 

along  this  astis,  and  iience  tlie  second  term  on  t>ie  right  in  I 

In  order  to  complete  the  proof,  we  need  only  snow  th-' 

is  bounded  as  a function  of  t and  z . S .jointly,  at  I'-.i.  t.  r 

small  t.  It  is  easy  to  show  tiiat  E ’J  (AZ^(t)  )' ] t,  E ' '^^t 

and  E*d  j ) ] /t  aj-e  boundeci  jointly  in  and  t.  TiiMs. 

Taylor's  Theorem  and  the  definition  of  , it  will  suffic'-  t ’ 


is 

T^'in  -1 
■ r ••.ft  i'-ieutl 
■lU  , t. 
from 

.:hov.'  that 


(f.^(/,)  E-AZj_(t)i  I-  f (r.)  E.^  Z>(t)l|.t 

is  bounded  jointly  in  t and  z.  Finally,  using  the  fact  tlia'. 


AZ<(tl  --  - va) 


and 


Z'-  't) 


X • 1 • Y • ( t 


for  t C 10,?“],  together  vit!i  the  botuuledness  of  the  secoii  j-oriler 
partial  derivatives  of  f ani  the  fact  that  (f,^  - f^l^v^.^.Oi  0, 

find  it  suffices  to  aiicw  that 


{ •.  E',[Y-(t)l] /t  (z.  E-d(M'(t)  - z ,)*!|/t 


Is  bounded  Jointly  in  t and  z = (z^,z).  hut  this  is  easy  to  verify 
from  the  known  first  passage  and  maximum  distributions  of  (one-dimensional. 
Brownian  Motion,  so  tiie  proof  is  complete. 


6.  A First  Passage  Problem 

We  now  show  that  l>(z)  - P'{T'  ■-  ■ I.  Thus,  h terminal--  ttie  Limit 

distribution  of  Z is  equivalent  to  solving  a first  ()asEa(:e  problen.  for 


Theorem  1.  “ P'[T’  '•  t)  for  t ^ 0 and 


Corollary.  t'(z  - l'^(T'  ■ for  z ..  S. 


Proof.  We  first  rewrite  Prop^osition  0 as 


(18)  ■■  Pq{Z;l  ■ -'l  '*■  - 


Next  obseT*ve  that 


- M,(t' 


sup  (X^(u)1  inf  ( , I - X^(ul I 


0-'  U'  t 


0 • u t 


.'0 


r 


(^0)  r - M,(t)  ^ ^ 


3UP  IX,,  (v)  ^ X (ul  I 
0 ■ V u • t 


inf  ( 1 
0 ^ V • u • L 


x,.(v)M  1-^^'^’-' 


Thus,  fief  ini  ;i<£ 


MKt;  - inf  (X'(ui]  , t _0  , 

^ 0 u t ^ 


M‘(t)  - inf  [X,'(v)  ^ X^(u)l  , t O . 

0 - V ■ u ' t 


v;e  comhine  (lS)-(;iO)  to  obtain 


(21)  4>^(z)  = 0,  M,*(t)  > 0]  - A M 


Finally,  ot'serve  that 


M^(t)  A M,*,(t;  inf  fX;(u  A X»(u)  »■  inf  X-'v  !l 
^ 0 - u ■ t ‘ ^ (IV. 


^ inf  IX' (u)  - Y*(uf|  • Inf  : 'j  ' . 

0 ^ u _ r,  ‘ 

Combining  this  with  (f'l)  and  the  definition  of  T'  eot:;;  tiie  jn-oof 

of  the  theorem.  The  corollary  follows  from  the  fact  Uiat  i ■^(  ) 

as  t t oo,  (See  Section  *0.  This  completes  the  pi'oof. 

Nov/  let  be  the  set  of  continuous  functions  f; .’>  * TR  whici. 

satisfy  the  followirif,  two  conditions. 

: i 


(22)  ail  partial  derivatives  of  f exist  and  are  continuous  evei7where 


on  S except  (possibly)  at  the  origin. 

(23)  Off  of  every  open  neighboi-hood  containing  tlie  origin,  f and  all 
of  its  parti<al  derivatives  are  bounded. 

From  (23'i  and  the  continuity  of  f it  follows  tint  f i.-  bounded. 

Theoi*em  Assume  0 and  + p,  ' C.  Sujdjoce  f C ^ satisfies 

(241  (f  - l(z, .0)  - 0 for  s,  0. 

4 • 1 i L 

(23)  D''?  = 0 excepit  (possibly  f.t  the  oriciii, 

(2(0  i'(7.^,z,  1 ~^0  unifomly  in  as  l •. 

('^7i  y(z^,z,)  ->  1 uniformly  in  ar  t <»  . 

Then  T = t. 

Remark  1 . There  is  rtroiii:'  reason  to  believe  tliat  the  conveise  is  al.so  true, 
meaning  that  'll  necessarily  satisfies  (:'4)-(  /')•  Fi’om  the  definition  ((1, 
it  follows  that  f satisfies  {-'(j]  aiii  (."Ij.  It  also  follows,  from  (1  ) 
that  ^ satisfies  (.4)  .so  long  as  both  partial  derivatives  exist.  Finally, 
from  the  corollary  to  Theorem  1 it  follows  easily  that  i>  and 

0*1  ■ 0.  To  prove  that  i ■;atisfieE  (.“  ) , it  thus  remains  only  to  show 
that  '1  belongs  to  a class  of  reasonably  smooth  functions  (like  f 


such  that  D'«f  is  bounded)  to  whicti  Proposition  ( ('d’t'  h'l')  can 

be  extended. 


Remark  It  is  important  to  note  that  Theorem  " requii-e.  • 

but  not  K <.  £ . In  Section  we  present  an  exajiiplf-  where  I-  is  in 
and  satisfies  (24' -{'-’Y)  but  where  all  partial  dex’ivatives  ol'  ; are 
unbounded  In  the  neighborhood  of  the  origin. 


Proof.  Suppose  0 < . <1.  and  assimie  i < ZJ'(O)  • lA  . Let  T be  the 
first  time  that  hits  either  , or  lA  . Then  T i ■ a stoppint’: 

time  for  Z" . and  it  is  easy  to  show  K^(T  ) >»  foi-  aid  such  that 

V < 7.J  < 1 ' . Proceeding  exactly  as  in  the  proof  of  Propu.sition  ' . one 
can  sho\-m  that  the  function  restricted  to  the  strip  ■ 1 

is  in  the  domain  of  tiie  generator  of  the  process  , inodifieii  by  absori  iiit 
barriers  at  ; and  l/>  . Moreover,  G'i'  - D'i  - (J,  where  i>"  now  denote; 
the  generator  of  the  modified  rrocerr.  Applying,  I'ynkin's  fonmila.  we 
then  have 

T 

(pH)  E"! ':'fz*(T  ))1  - 'i'(z)  + E-  / G>r(ZM't))  dt  ;'(■,)  . 

; . 2 0 

Observe  that  T ->T‘  almost  sui-ely  as  ■ i 0.  Also,  frotu  ' . anu 

(,•7)  it  follows  tiiat  t.he  Left  side  of  (pH)  converges  to  i -iT'  : w) 
as  t 0.  Thus,  letting  1 0 in  (:’H)  anci  usinj’;  ti.e  Corollary  to 
Theorem  1,  we  obtain  ^(s,)  'S‘(zi, 


7.  T^ie  Density  l■^InGtlon 

The  remainder  of  the  paper  i-onaii^ts  of  detuilei  ou ' alntior...  for 
which  subscripted  variables  are  clumsy.  Thus  {)oint;-  in  will  be  ienor’:! 
(x,y)  hereafter,  rather  that  ^ . • Let  b-  tf.e  set 

of  functions  f;S  -»  Di  that  are  continuous  except  fpossihly)  et  the  oricin 
and  that  .-atisfy  {c’j)  and  ( '5).  A function  in  ^ may  be  unbour.  leu  li; 

the  neighborhood  of  the  origin. 

Theorem  Assuise  < 0 and  t-  < 0,  Tliere  exi.sts  a 1-unction 

f r Uq  satisfying  iff  there  exists  f ^ ’ satisfying 

(P9)  D'-f  - 0 on  the  int.erior  of  G . 


(50) 

t i a-'f  . - p fl(x.O)  ■ 0 

for  x 

(51) 

'1  ‘^l^’l  ‘’l  ' 'd  ■ ■ l'"" 

. V t 0 

for  y 

(5^) 

JU  JU 

f > 0 and 

f(x.y  ‘ ix  dy 

1 . 

0 0 

In  this  case,  f = - i'  , and 


X X'l  y-u 

(5^)  'V{z,y)  f f f(u,vi  dvdu  for  (x,y'  t G. 

0 0 


Corollary.  If  f .-.atisfies  ( then  f is  a density  for  the 


'U 


llinit  distribution  F. 


Remark  1.  The  comments  following  Tlieorem  when  combined  v/iU.  Theorem 
suggest  that  (;.•"')-(■■.  ) Hva  h(;Ui  u«oe.'.L.ai^  ami  aaffi.-ic-ni  i'~ir  Uic  dunoity  of 
the  Limit  distribution. 


Remark  ih  Again  v/e  emphaoi  .;e  that  the  theorem  doeo  not  regui-r-  f 
In  the  example  of  fieetion  ■,  both  f and  its  pa'-ti  il  ..ori  ver, 
are  unboun  ie*i  iii  the  neif'hborfiood  of  ti.c  origin. 


Proof.  First  suppo.se  that  'if  satisfies  ' - ' aiH  lefine 

f ‘ excepc  at  the  origin  (where  the  value  of  f ui  ije  .,et 

arbitrarily).  From  the  definition  of  it  foUowr.  iiiai  f , 

Also,  since  (T.^  - t )(x,0)  0,  we  have 


^ 2-X 

(3*t^  v^'J'(x,?.-x)  (i'^  - -I'.KxjZ-xi  ---■  ; fvx.v)  dv 

0 


for  z > A ' 0.  rje.xt,  since  v(O.z)  0,  we  liave 


(,5b) 


X ^ X ;;-n 

'i’(x,z-x','  " - — ?fu,:--u'  du 

^ 'U 

<1  0 0 


for  z > X 0,  and  is  enuivalen*.  to  ■ . N ■ fr..in  ; ,,'e  directly 

compute  that  p-  ' f on  tne  in',  rior  of  d.  wnere 


1 

Vfi'x.yi  D:f(y.v'(iv  i r-f(u.yhlu  ! {-  ei,  ^ a,  ' f(x.yi. 


i " 1 


Djf  “ interior  of  fi, 


(i  0-'^  t-  0'^  ■^  h o'.'f':,  ' ^'1  ■*■ 


oil  the  interior  of 


From  the  fact  that  D'i'  - 0 on  the  interior  of  0 v/e  have 

.1.'  ,2  ^2 

(56)  0 - ~)  D^Kx.y)  = - -—)  .^f(x,y  , 

■y  'y 


(57)  0 = (^  - ^)  i>-Kx,u)  - C 


) A"  i'  ( X , 0 ) , 


(58)  0 D^'i(0,y)  - An'(0,y), 


Directly  computing  the  right-most  expressions  in  (56)- '5',.  fro:i.  (55  . a-e 
obtain  (2ti)-(51).  Finally,  Theorem  • gives  Y - 5,  and  then  frimi  ' '■ 
and  the  definition  (h)  of  'i  it  follows  that  f is  a density  for  the 
proper  distribution  F.  Thus  f satisfies  i,52). 

For  the  converse,  suppose  f satisfies  an.:  lefiue 

ty  (55).  It  is  immediate  from  (52)  and  (55)  tliat  ■'  sati-d'ie;:  ! , 

(;.'6)  and  Thus  \re  need  only  show  D*T  0.  Again  we  nav-  t ■■■  -('  ' 

the  second  equality  in  each  case  following  from  the  fact  tluil  ' ' f 

by  (55),  and  the  triple  of  first  equalities  being  eouivalent  to  (.  '*  -^It. 
But  (5<  )-('"8'.  ai'e  necessai’y  and  sufficient  for  h‘ • K (a  con.  tanl  . 
Thus  we  need  to  shovi  K - 0.  Wiiatever  the  value  of  K,  t!ie  proof  of 

Theorem  , need  be  all.ereJ  only  trivially  to  establi.Ui  ttial. 


E;',  *(T  ))l  '>(..)  ' KE;,iT  . 


Wlien 

< 0 and  1 

, <0,  it  is 

easy 

to  show 

that  E'(T  ) O'. 

as  . 

i 0 whenever 

^ 0.  Oince 

i'‘  is 

bou  tided, 

it  i'ulluws  from  ( '/ 

tliat 

K r.  0. 

H.  Exponential  Solutions 

lUnoe  our  banic  i if ferential  equation  D'f  0 Im.;  courtant 
coefficients,  it  is  natural  to  seek  solutions  of  t)je  se[.a/-abie  exi.onential 
form 

(»'0)  f(x.yi  . :x.y)  ( ,i. 

Proposition  !.  Assume  0 and  ■-  Cu  'riiere  exi.-.t,.-  a solution 

f of  (p'di-fiPl  having  tiie  fom  (>t0'<  iff  ~ + Sj  in  liioh  case 

<v  - Pluii  /Pj  and  (-  - V . 

Remark.  A solution  f of  ( can  have  tin?  fnnii  f ;-;.yi  k(x‘  Iv'y 

only  if  f is  of  tlie  forrr,  (-'tO'.  This  follows  easil;,  from  t,;io  boundary 

conditions  (iO)  and  fill.  Thus  the  proposition  implies  ins'  -o'’  i i.  ' 0 

L L 

is  necessary  anti  sufficient  for  the  existence  of  a .'eparalue  .-oliition. 

Proof.  Assuming  the  form  fliO),  conditions  (P'fi-f-,1  are  equivalent  to 
(^<1  ) ^ ^ ~ e h'  + ,h  <">  . 


k + hp  = 0 ’ 


\ '^V  ' ^ °'l  * ‘^'l'  *■  M]_  = 0 


Multiplying!;  (42^  througii  by  ('rfhich  must  be  positive  if  f is  to  be  a 
density)  and  subtracting  this  from  (4l).  we  obtain  - o'^''  + Comparing 

this  with  (4^),  we  see  that  (4l)-(4ji)  can  hold  simultaneously  iff 
^ CTf  + cJ^,,  0,  in  which  case  u - : Jub.-titutinf.:  this  value 

of  a in  (4?)  gives  , 2]^^  + u^l/a',. 


9.  Another  Special  Case 

Wiien  the  density  f does  not  have  the  simple  form  computirv: 

it  explicitly  appears  to  be  a difficult  problem,  Oply  one  other  speci!-! 
case  has  been  solved  thus  far. 


Proposition  ii.  Assume  cr'^,,  = 0,  o'^  ■ a,]  - 1,  pp  ^ 0,  and  •-  0.  Let 
M ~ IutI-  Then  the  .lensity  of  the  limit  distribution  K is 


f(x,y)  = Ke~^^l  r(x,y)  e cos  ~d(x.y  | , 


where  r(x,y)  (x'‘  + y . d(x,y)  - tan"^(y  x'  and  K 'r  ^ 


Remark . The  si.ightly  more  general  case  cr'^  = o',  o can  of  course  be 


accommodated  by  a rescaling. 


1^ 


Proof.  In  this  case  our  equations  (iry)  through  for  f become 

(hJ()  i f t-  i f I-  ,;f  - 0 on  the  interior, 

(45)  i f,  (x,0)  - 0 for  X ■'  0 

(4(  ) f.|  I i r.  I fi(0,y)  0 for  y . 0 . 

Letting  g(x,yi  - f(>:,.yi,  we  re-express  (44)  tiirough  (4i  in  terms 

of  g as 

(4?)  g !■  , - pg  on  the  interior, 

i.  L 

(4rt)  g,/x,0)  0 for  X 0 , 


Hypothesi<:ing  a solution,  of  the  form  h(r,ri)  = 9(r)  ^{O),  we  find  that 
_ 1 A) 

cp(r^  ^ r exp(-^.r)  and  - cos(0/:’)  satisfy  (v0)-(‘:-).  Kever.siry: 

the  transfoi-mations  then  gives  the  desired  formula  for  f except  that 
the  normalisation  constant  K has  not  been  deten.iined . We  need 


iff  f(x.y  <ixiy  K / / lr(x,y)]”^  ' e'*-  'cosl(l  '(x,: 


“ -1  • -,J  l-i-cos  ->lr  A 

- K / ' E e ' Gos(—  i9)r  Jrclt?  . 

0 0 


Let  X(t?)  ^ ^i(  1 i-  cor,  OJ,  go  X(0)  - -',i  cos  (—  #9).  Observe  tiiat 


f ir  . /’  e'"  du  ^ r(^) 


~ vf  A"''  “(3)  - ~ n ‘ cos  (--))  r 


^ V « ( cog"^(7  ;)) 


Thus  we  need 


rri  lA/o  T 


1 K[i  / cos''’(i  d)d.9  - K.i 


i l/:\ 


tari(-  J 


1 


10 , Concludlnfi  Remarkr 

Probably  the  best  known  work  on  mul  ti -dimensional  d i ‘‘riK-lonr-;  with 
boundaries  is  the  seminal  paper  by  Stroock  and  Varadluin  ■>  !■  Otriotiy 

speaking,  our  process  Z lies  outside  the  Stroock-Varadhan  theci-y,  because  j 

its  boundary  fails  to  meet  their  smoothness  requirements.  Tne  sanie  in  true 

of  the  theory  developed  by  Watanabe  1 ('  1 • | 

I 

i 

In  Section  1 we  have  discussed  only  the  fact  that  the  equilibrium  ’ 

distribution  of  the  diffusion  Z approximates  the  disti'lbution  of  j 

! 

for  a tandem  uue\iint:  system,  where  W is  the  equilibrium  'waiting  time  | 

vector.  One  can,  ho-wever,  shov;  that  in  heavy  traffic  the  li  sti-ibution  of  j 

the  entire  (vectorl  waiting  time  process  is  approximated  Ly  the  distribution  | 

of  Z in  the  following',  sense.  Assume  a stable  tandem  (;uenii;g  system 
(u  >0),  and  define  a two-dimensional  process 

f (t  ; = , t > 0 . I 

where  1 • ] denotes  integer  paru.  Now  consider  a seauerice  of  stable  j 

tandem  systems  with  S'  i 0.  Under  a.ssuinptions  like  those  emj  loyed  in 

our  previous  paper  [wj,  and  using  very  tiaich  the  same  methods,  it  can  be  ’ 

shown  that  ' conveiges  weakly  (in  the  appropriate  function  sp.ace)  to 
Z as  •'(  t 0.  In  a similar  vein,  one  mignt  examine  the  nonnaii.ed 
(vector)  queue  length  process 

iy(t)  ■-  <<Q(t,'  j""! , t > 0 . 

•where  Q(ti  Qj^(t),  Q (t)  ] an.i  denotes  the  numoer  of  customer^ 
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present  at  station  k at  time  t.  Igletiart  and  Wiiitt  I i.ave  pj  oved 


heavy  traffic  limit  theoreii.;;  foi'  processes  of  tliis  I'onn  as^oelated  v/itii 


general  class  of  acyclic  queuing  networks.  (Tlieii  nonnal  j at  ions , 


however,  are  expressed  in  a rattier  different  fom.  ) If  one  speciali.'.et 


their  results  for  i-  to  the  case  of  two  queues  in  tandem,  one  finds 


that  converges  weakly  to  Z as  '<  1 0,  where  Z is  e„.  eni.ial iy  the 


same  difl>ision  studied  here. 


The  last  paragrap-h  .suggests  that  we  might  be  interesLe  i in  quantities 


associated  vtith  Z other  than  its  equilibri'um  distribution.  In  particular. 


the  quantity  provides  an  approx limti on  for  tne  tran  ^ lent  waiting 


time  or  queue  leni_ith  (iistribution  of  a tandem  system,  provi  iei  that  the 


system  stai’ts  empty  and  heavy  traffic  conditions  prevail,  iiieorem  1 


and  Proposition  1 together  suggest  tliat  should  sati.  t'y 


c)  - > 0 i-ii  i-ho  irit(.‘rioi'  of  .S, 


together  with  the  boundary  conditions  'l'.(0,z,,)  ~ 0 and 


't>,  .0)  =0  foi  t ■■  0 and  , ■ 0 . 

oz,  t'  1 1 


.and  the  initial  con.iition  '!’q(')  = 1-  We  .sliall  iiuil.e  no  attemi  t to  pr.ve 


thl.s,  except  for  the  folJ.ov'ing  comment.  If  a .sufficiently  .pootli  .solution 


can  be  found,  then  uynkin's  formula  can  be  used  to  prov>.-  tlial  ' bis  .lolution 


is  f • 1 , .just  as  in 


the  proof  of  llieorem 


I 


m 
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ABSTRACT 


, Consider  a.  pair  of  single  server  (jueues  arranged  In  series.-  (Tiiis 

j is  the  simplest  example  of  a queuing  netv(ork.  ) in  an  earlier  paper  [2], 

—i  a limit  theorem  vras  proved  to  justify  a heavy  traffic  approximation  for 
I the  (two-dimensional)  equilibrium  v/aiting  time  distribution.  Specifically 
the  v^aiting  time  distribution  -uas  shovm  to  be  approximiated  by  the  limit 
^ distribution  P'  of  a certain  vector  stochastic  process  Z.  Tlie  process  Z 
•rfas  definea  as  an  explicit,  but  relatively  complicated,  transfoimation  of 
vector  Brownian  Motion,  and  the  general  problem  of  detern;inir,g  F was 
left  unsolved. 


I 

i 

I 


It  is  shonn  i-rr-ahis  paper;,  that  Z is  a diffusion  process  (con- 
tinuous strong  Markov  process)  whose  state  space  S is  the  non-negative 
qua.lrant.  On  the  interior  of  S,  the  process  behaves  as  an  ordinary 
vector  Brownian  Motion,  and  it  reflects  instantaneously  at  each  boundary 
surface  (axis).  At  one  axis,  the  reflection  is  normal,  but  at  the  other 
axis  "it  has  a tangential  com.ponent  as  well.  Tlie  generator  of  7,  is  cal- 
culated. 


It  is  shown  tliat  the  limit  distribution  F is  the  solution  of  a 
first  passage  problem  for  a certain  dual  diffusion  process  Z*-.  The  generator 
of  Z'  is  calculated,  and  the  analytica''  theory  of  Markov  process  is  used 
to  lei'lve  a partial  differential  equation  (with  boundary  conditions)  for 
the  density  f of  F.  Necessary  and  sufficient  conditions  are  found  for 
f to  be  separable  (for  the  Limit  distribution  to  have  independent  com- 
ponents).^ This  exteiid.-,  slightly  the  class  of  explicit  solutions  found 
..lovlousiy  vr,  ’ ;.  Another  special  case  is  solved  explicitly,  showing 

that  the  densi'ty  i.s  not  in  general  separable. 
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